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Ground-state properties of 16 O and 40 Ca are calculated with a low- momentum nucleon-nucleon 
potential, VJow-kj derived from the chiral N 3 LO interaction recently constructed by Entem and 
Machleidt. The smooth Vi ow -k is used directly in a Hartree-Fock approach, avoiding the difficulties 
of the Brueckner-Hartree-Fock procedure. Corrections up to third order in the Goldstone expansion 
are evaluated, leading to results that are in very good agreement with experiment. Convergence 
properties of the expansion are examined. 
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I. INTRODUCTION 

A fundamental problem in nuclear theory has long 
been the calculation of the bulk properties of closed-shell 
nuclei, such as their binding energy and charge radius, 
starting from realistic nucleon-nucleon (NN) potentials. 
Potentials like CD-Bonn p|, Nijmegen 0], Argonne vi$ 
@ , and the new chiral potential of Ref. (J reproduce the 
NN scattering data and the observed deuteron proper- 
ties very accurately, but, because of their strong repulsion 
at short distances, none of them can be used directly in 
nuclear structure calculations. 

A traditional approach to this problem is the 
Brueckner-Goldstone (BG) theory where the Gold- 
stone perturbative expansion is re-ordered summing to all 
orders a selected class of diagrams, the ladder diagrams. 
This implies replacing the bare interaction (Vat at) ver- 
tices by the reaction matrix (G) ones and omitting the 
ladder diagrams. Within this framework, one has the 
well known Brueckner-Hartree-Fock (BHF) theory when 
the self-consistent definition is adopted for the single- 
particle (SP) auxiliary potential and only the first-order 
contribution in the BG expansion is taken into account. 
The BHF approximation gives therefore a mean field de- 
scription of the ground state of nuclei in terms of the 
G matrix, the latter taking into account the correlations 
between pairs of nucleons. However, owing to the energy 
dependence of G, this procedure is not without difficul- 
ties. An important issue is the choice of the SP energies 
for states above the Fermi surface, which are not uniquely 
defined @. 

Calculations for finite nuclei within the BHF approach 
lead usually to insufficient binding energy as well as too 
small charge radii []]]. In this context, extensions of the 
conventional BHF approach have been proposed to ac- 
count for long-range correlations Q • It is worth mention- 
ing that alternative approaches have been developed, as 
for instance the correlated basis function method y| and 
the coupled cluster method ^} > m both of which corre- 



lations are directly embedded into the wave functions. A 
comprehensive review of the various methods is given in 
Ref. ' which also includes a discussion of calculations 
for nuclear matter as well as for finite nuclei. An other 
effort in this direction is the uni tary model-operator ap- 
proach of Suzuki and Okamoto ■ 

Recently, a new technique to renormalize the short- 
range repulsion of a realistic NN potential by integrat- 
ing out its high momentum components has been pro- 
posed [H[l|. The resulting low-momentum potential, 
which we call Viow-k, is a smooth potential that preserves 
the low-energy physics of Vnn and is therefore suitable 
for being used directly in nuclear structure calculations. 
We have employed 14 ow -k derived from modern NN po- 
tentials to calculate shell-model effective interactions by 
means of the Q-box plus folded diagram method. Within 
this framework, several nuclei with few valence particles 
have been studied 0, 0, , leading to the conclusion 
that Mow-k is a valid input for realistic shell-model cal- 
culations. 

The main purpose of this paper is to show that this 
potential may also be profitably used for the calculation 
of ground state properties of doubly closed nuclei. To this 
end, we have employed the Goldstone expansion which, 
given the smooth behavior of Viow-k, does not require any 
rearrangement. As a first step of our procedure, we solve 
the Hartree-Fock (HF) equations for Vj ow -k- Then, using 
the obtained self-consistent field as auxiliary potential, 
we calculate the Goldstone expansion including diagrams 
up to third order in Viow-k ■ 

Here, we present the results obtained for 16 O and 40 Ca 
starting from the new realistic NN potential of Entem 
and Machleidt Q based on chiral perturbation theory at 
the next-to-next-to-next-to-leading order (N 3 LO; fourth 
order). This potential is an improved version of the ear- 
lier chiral NN potential (known as Idaho potential) 
constructed by the same authors, which includes two- 
pion exchange contributions only up to chiral order three. 
We have recently employed the latter potential in shell- 
model calculations for various two-particle valence nuclei 
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The paper is organized as follows. In Sec. II we first 
describe the main features of the derivation of Mow-k, 
then outline the essentials of our calculation. In Sec. 
Ill we present our results and compare them with the 
experimental data. Some concluding remarks are given 
in Sec. IV. 



II. METHOD OF CALCULATION 

The first step in our approach is to integrate out 
the high- momentum components of Vnn- According to 
the general definition of a renormalization group trans- 
formation, the decimation must be such that the low- 
energy observables calculated in the full theory are pre- 
served exactly by the effective theory. Once the relevant 
low-energy modes are identified, all remaining modes or 
states have to be integrated out. 

For the nucleon-nucleon problem in vacuum, we re- 
quire that the deuteron binding energy, low-energy phase 
shifts, and low-momentum half-on-shell T matrix calcu- 
lated from Vnn must be reproduced by Viow-k- 

The full-space Schrodinger equation may be written as 



= H = Hq + Vnn, 



(1) 



where Hq is the unperturbed Hamiltonian, namely, the 
kinetic energy. The above equation can be reduced to a 
model-space one of the form 

PHckP^h = £ M P1V; H eS = H Q + Mow-k, (2) 

where P denotes the model-space, which is defined by 
momentum fc < fc cut = A, fc being the relative momentum 
and fc cut a cut-off momentum. 

The half-on-shell T matrix of Vn n is 



T(k',k,k 2 ) = V NN (k',k) + 



q 2 dqV N N(k',q)- 



1 



-T(q,k,k 2 ) , (3) 



/ k 2 -q 2 + iO+~ 

and the effective low-momentum T matrix is defined by 



[l3L fill it has been shown that the above requirements 
are satisfied when Viow-k is given by the folded-diagram 
series 



= Q-Q' / Q 



-Q' Q Q Q+ 



Mow-k 

(5) 

where Q is an irreducible vertex function, in the sense 
that its intermediate states must be outside the model 
space P. The integral sign represents a generalized fold- 
ing operation |19| , and Q' is obtained from Q by removing 
terms of first order in the interaction. 

The above V[ ow _k can be calculated by means of itera- 
tive techniques. We have used here an iteration method 
proposed in Ref. |2C|. which is particularly suitable for 
non-degenerate model spaces. This method, which we 
refer to as Andreozzi-Lee-Suzuki (ALS) method, is an it- 
erative method of the Lee-Suzuki type, which converges 
to the lowest d eigenvalues of H 1 d being the dimension 
of the P space. Since the Vj ow -k obtained by this tech- 
nique is non-hermitian, we have made use of the simple 
and numerically convenient hermitization procedure sug- 
gested in Ref. [2(|. We have verified that the deuteron 
binding energy and the phase shifts up to the cut-off mo- 
mentum A are preserved by Vjow-k- 

An important question in our approach is what value 
one should use for the cut-off momentum. A discussion 
of this point as well as a criterion for the choice of A can 
be found in Ref. According to this criterion, we 

have used here A = 2.1 frn~ . 
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Ti ow ^(p',p,p 2 ) = V ow - k (p,p) + 



/" 

Jo 



q 2 dqV ow - k (p',q)- 



P 



i0+ 



T low _ k (g,p,p 2 ) . (4) 



Note that for Xi ow _k the intermediate states are inte- 
grated up to A. 

It is required that, for p and p' both belonging to 
P (p,p' < A), T(p',p,p 2 ) = Ti ow _ k (p',p,p 2 ). In Refs. 



FIG. 1: Behavior of E KF with huj and N for lfa O. 

We have found it appropriate, however, to check on 
the sensitivity of our results to moderate changes in the 
value of A. It has turned out that they change very little 
when letting A vary from 2.0 to 2.2 fm _1 . For instance, 
at the second order in the Goldstone expansion (see next 
section), the binding energy per nucleon for 16 is 7.29 
and 7.12 MeV for A = 2.0 and A = 2.2 fm _1 , respectively. 

As already mentioned in the Introduction, our starting 
point is the N 3 LO potential |4|. From this potential we 
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derive the corresponding Mow-k and use it directly in a 
HF calculation. The HF equations are then solved for 
16 O and 40 Ca making use of a harmonic-oscillator basis. 
We assume spherical symmetry, which implies that the 
HF SP states \a) have good orbital and total angular 
momentum. Therefore, they can be expanded in terms 
of oscillator wave functions 
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FIG. 2: Behavior of £ H p with hw and TV for 40 Ca. 



2, where we show the behavior of the HF ground-state 
energy of 16 O and 40 Ca versus hw for different values of 
N. The results for N — 5 are quite stable. In our cal- 
culations the values of ftw have been derived from the 
expression fmj = 45A~ 1 / 3 — 25 A~ 2 I Z [2l|. which repro- 
duces the rms radii in an independent-particle approx- 
imation with harmonic-oscillator wave functions. This 
expression gives Tilo = 14 and 11 MeV for 16 and 40 Ca, 
respectively. In solving Eq. (7), the Coulomb potential 
has been added to Vj ow _k for protons. 




(c) (d) (e) 



FIG. 3: First-, second-, and third-order diagrams in the Gold- 
stone expansion. 



(6) 



where the sum is over the principal quantum number 
only. The expansion coefficients are determined by 
solving self-consistently the HF equations 



J2W + W)C%=e a CZ 
a' 



(7) 



where t is the kinetic energy and the HF potential U is 
defined as 



{n\U\n') = y^^a h \Vi ow ^ k \^a h ) 



(8) 



with the index referring to occupied states in the 
ground-state HF Slater determinant. 

Once Eq. (7) has been solved, the ground-state prop- 
erties of the nucleus can be calculated. In particular, the 
total energy has the well-known expression 



-^hf = \ [( a h\A a h) + e Q J 



(9) 



In our calculations the sum in the expansion (6) has 
been extended up to N = 5 terms. We have verified that 
this truncation is sufficient to ensure that the HF results 
do not significantly depend on the variation of the os- 
cillator constant hui. This is illustrated in Figs. 1 and 



We use the HF basis to sum the Goldstone expansion 
including contributions up to third order in Viow-k- In 
Fig. 3 we report the first-, second-, and third-order 
diagrams |22|. The intermediate states involved in the 
evaluation of these diagrams are those obtained from the 
lowest 10 and 11 h.o. major shells for 16 and 40 Ca, re- 
spectively. Our results show convergence for these large 
spaces: for example, the 16 binding energy per nucleon 
in second order approximation is 7.03, 7.22, and 7.30 
MeV when considering 9, 10, and 11 major shells, re- 
spectively. Similar diagrams have been used to calculate 
the corrections to the rms radius. As is well known, re- 
taining only the first-order term in this expansion yields 
the HF results. 

To conclude, it is worth stressing that using the Gold- 
stone expansion in terms of Haw-k one has to include 
also diagrams (b) and (c). This is not the case of the BG 
approach, where these diagrams are already contained in 
diagram (a) through the G matrix. 



TABLE I: Comparison of the calculated binding energy per 
nucleon (MeV/nucleon) and rms charge radius (fm) with the 
experimental data for 16 O and 40 Ca. 



Nucleus 



HF HF+2nd HF+2nd+3rd 



Expt. 



'O 



J Ca 



B/A 3.23 
(r e ) 2.30 



7.22 
2.52 



B/A 6.19 9.10 
(r c ) 2.610 3.302 



7.52 
2.65 

9.19 
3.444 



7.98 
2.73 ± 0.02 

8.55 
3.485 ± 0.003 
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III. RESULTS 



In Table I, the calculated binding energy per nucleon 
and the rms charge radius for both 16 and 40 Ca are 
compared with the experimental data [2^, I24I |25| . This 
Table contains the HF results as well as the values ob- 
tained including second- and third-order contributions. 
We see that the renormalization of the short range re- 
pulsion through Viow-k is sufficient to yield positive HF 
binding energies, albeit too small as compared to the 
experimental values. We also see that the HF approxi- 
mation significantly underestimates the rms radii. This 
evidences the role of higher-order contributions in the 
Goldstone expansion to account for correlations beyond 
the mean field. The binding energies and radii calculated 
including diagrams up to third order are very satisfac- 
tory. In fact, the binding energies gain about 4 and 3 
MeV for 16 O and 40 Ca, respectively, while the radii in- 
crease by about 0.4 and 0.8 fm coming quite close to the 
experimental values. 



TABLE III: Calculated and experimental single-particle en- 
ergies (MeV) for 40 Ca. Experimental data are taken from 
[2ll2ll27t. 

Neutron Proton 



Orbital 


Calc. 


Expt. 


Calc. 


Expt. 


Sl/2 


-97.944 




-87.501 


-49.1 ± 12 
-77 ± 14 


PS/2 


-63.760 




-53.934 


-33.3 ±6.5 


Pl/2 


-54.959 




-45.269 


-32 ±4 


d*,/2 


-33.018 


-21.30 


-23.749 


-14.9 ±2.5 
-13.8 ±7.5 


Sl/2 


-27.406 


-18.104 


-18.238 


-10.850 


^3/2 


-19.595 


-15.641 


-10.663 


-8.328 


/t/2 


-6.579 


-8.363 


2.047 


-1.085 


PS/2 


-4.325 


-6.420 


3.293 


0.631 


Pl/2 


-0.973 




5.865 




h/2 


5.852 




12.484 





IV. CONCLUDING REMARKS 



TABLE II: Calculated and experimental single-particle en- 
ergies (Me V) for le O. Experimental data are taken from 







16Q 








Neutron 


Proton 


Orbital 


Calc. 


Expt. 


Calc. 


Expt. 


Sl/2 


-53.786 


-47 


-48.606 


-44 ± 7 


P3/2 


-23.225 


-21.839 


-19.264 


-18.451 


Pl/2 


-15.649 


-15.663 


-11.841 


-12.127 


ffc/2 


-0.056 


-4.144 


3.564 


-0.601 


Sl/2 


-0.481 


-3.273 


2.651 


-0.106 


^3/2 


5.814 


0.941 


8.713 


4.399 



A discussion of the convergence properties of the per- 
turbative series is now in order. To this end, the HF 
potential energy, the second- , and third-order correc- 
tions have to be compared. In 16 the HF potential 
energy per nucleon is = —23.5 MeV, which is ob- 
tained by subtracting from the total HF energy the con- 
tribution of the kinetic term. Thus for the ratio of the 
second- to first-order term we obtain V {2) /V^ = 0.17, 
while the ratio /V™ is 0.08. Similarly for 40 Ca 
we have V^> = -33.7 MeV, V™ /yU) = 0.09, and 
v (3)/y{2) = 03 on these grounds, we may conclude 
that the convergence of the series is fairly rapid and that 
higher-order contributions are negligible. 

For the sake of completeness, in Tables II and III we 
report the HF single-hole energies as well as the energies 
of the low- lying particle states of 16 O and 40 Ca. In Ta- 
bles IV and V the calculated occupation probabilities for 
states up to the Fermi level are reported. 



In this work, starting from the chiral N 3 LO poten- 
tial 0, we have performed calculations for the ground- 
state properties of the doubly closed nuclei 16 O and 40 Ca 
making use of the Goldstone expansion. This has been 
done within the framework of a new approach [l3l Il4| 
to the renormalization of the short-range repulsion of re- 
alistic NN potentials, wherein a low-momentum poten- 
tial Viow-k is constructed which preserves the low-energy 
physics of the original potential. We consider a main 
achievement of our study to have shown that How— k is 
suitable for being used directly in the Goldstone expan- 
sion. Namely, unlike the traditional BHF approach, there 
is no need to first calculate the G matrix. We have seen 
that taking into account higher-order contributions (es- 
sentially the second-order terms) of Vjow-k in the Gold- 
stone expansion yields very good results for the binding 
energy and charge radius of 16 and 40 Ca. 



TABLE IV: 


Calculated occupation 


probabilities for le O. 


Orbital 


Proton 


Neutron 


Sl/2 


0.881 


0.881 


P3/2 


0.822 


0.818 


Pl/2 


0.767 


0.760 



In this context, it is worth emphasizing that the idea 
of bypassing the G-matrix approach to nuclear structure 
calculations is not new [2^, [23. From this viewpoint, the 
present study comes close in spirit to the early work car- 
ried out at the MIT HHHEll in the mid 1960's. There 
the Tabakin's separable nonlocal NN potential [2^] was 
used in HF calculations and then the second-order correc- 
tion to the binding energy was evaluated, obtaining quite 
satisfactory results. It is indeed very gratifying that our 
results based on the use of the Vjow-k derived from the 
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modern N 3 LO potential come significantly closer to the 
experimental data. It should be noted that our results are 
quite good also when compared to those of recent BHF 
calculations 8], where different modern NN potentials 
have been used and long-range cortrelations have been 
considered within the framework of the Green function 
approach. 



In summary, we may conclude that the results of the 
present study together with those of our recent shell- 
model calculations, show that the V[ ow -k approach pro- 
vides a simple and reliable way of "smoothing out" the 
repulsive core contained in the modern NN potentials 
before using them in microscopic nuclear structure cal- 
culations. 



TABLE V: Calculated occupation probabilities for c Ca. 



Orbital Proton Neutron 

~~ sITa 51)45 0.947 

p 3/2 0.931 0.932 

Pi/a 0.921 0.922 

d 5/2 0.885 0.884 

si/2 0.858 0.855 

d 3/2 0.811 0.807 
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